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We present an analytical study of a nonlinear oscillator subject to an additive Ornstein-Uhlenbeck 
noise. Known results are mainly perturbative and are restricted to the large dissipation limit 
(obtained by neglecting the inertial term) or to a quasi- white noise (i.e., a noise with vanishingly 
small correlation time). Here, in contrast, we study the small dissipation case (we retain the inertial 
term) and consider a noise with finite correlation time. Our analysis is non perturbative and based 
on a recursive adiabatic elimination scheme: a reduced effective Langevin dynamics for the slow 
action variable is obtained after averaging out the fast angular variable. In the conservative case, 
we show that the physical observables grow algebraically with time and calculate the associated 
anomalous scaling exponents and generalized diffusion constants. In the case of small dissipation, 
we derive an analytic expression of the stationary Probability Distribution Function (P.D.F.) which 
differs from the canonical Boltzmann-Gibbs distribution. Our results are in excellent agreement 
with numerical simulations. 
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I. INTRODUCTION 

A theoretical description of the intrinsic thermal fluctuations induced by the coupling of a given dynamical system 
with a heat bath can be obtained by adding a noise term to the equations of motion. From this point of view that 
originates from Langevin's treatment of Brownian motion, the effect of the heat bath is modeled as a random force 
that transforms a deterministic motion into a stochastic dynamics in phase space. The physics of the system is then 
described by a probability distribution function (P.D.F.) and the information about the observables is of statistical 
nature. When the additive noise is white (i.e., the correlation time of the noise is smaller than any intrinsic time 
scale of the system) and its amplitude satisfies the fluctuation-dissipation relation, the stationary P.D.F. is identical 
to the canonical Boltzmann-Gibbs distribution [HQ 

A paradigm of a stochastic dynamical system for a quantitative study of the interplay between randomness and 
nonlinearity is provided by a particle trapped in a nonlinear confining potential and subject to a random noise 
We showed in 0, Q that in the limit of a vanishingly small damping rate, such a particle exhibits anomalous 
diffusion with exponents that depend on the form of the confining potential at infinity. For a non-zero damping rate, 
this anomalous diffusion occurs as intermediate time asymptotics: the particle diffuses in phase space and absorbs 
energy from the noise until the dissipative time scale is reached and the physical observables become stationary. For 
an additive Gaussian white noise, it was shown explicitly Q that, for times larger than the inverse damping rate, 
the intermediate time P.D.F. matches the canonical distribution. For colored noise, we observed from numerical 
simulations that the anomalous growth exponents are different from those calculated for white noise. This behavior 
strongly suggests that the long time asymptotics in the case of colored noise is not identical to the canonical Boltzmann- 
Gibbs distribution. The physical reason for this change of behavior is the following: the intrinsic period of a nonlinear 
oscillator is a decreasing function of its energy and when the amplitude of the oscillations grows the period eventually 
becomes shorter than the correlation time of the noise. In that case, destructive interference between the fast variable 
and the noise suppresses the energy transfer from the noise to the system and the diffusion slows down. In this regime, 
the correlation time of the noise ceases to be the shortest time scale in the system and the noise can not be treated 
perturbatively as white. Therefore, usual perturbative calculations based on small correlation time expansions [sLIfillOjl 
cannot predict the correct colored noise scalings (this was shown explicitely for the pendulum with multiplicative noise 
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in Besides, the averaging procedure used for the white noise fails for colored noise and in [f| we could only infer 

the colored noise exponents using qualitative scaling arguments. 

Our aim, in this work, is to perform an analytical study of nonlinear oscillator subject to an additive Ornstein- 
Uhlcnbeck noise and to derive quantitative results for its long time behavior. We need to generalize the averaging 
method used for white noise in [a, [(J so that it can be applied to a regime where the correlation time of the noise is not 
the shortest time scale in the system. We define, in a recursive manner, new coordinates on phase space that embody, 
order by order, correlations between the noise and the fast angular variable and then average out the fast variable. 
This technique allows us to obtain an analytical expression for the long time asymptotic behavior of the P.D.F. of 
the system, to derive the anomalous scaling exponents and to calculate generalized diffusion constants. In particular, 
it is shown explicitly that the long time P.D.F. of a nonlinear oscillator subject to additive colored noise is not the 
canonical Boltzmann-Gibbs distribution. Our results are systematically compared with numerical simulations. 

The outline of this paper is as follows. In section [HJ we review results derived in and which are required 
for the present work. In section IIIII we first explain why the straightforward averaging technique fails for colored 
noise, provide a heuristic calculation that explains the idea underlying the correct averaging scheme and develop it 
analytically. An effective averaged Langevin system for the slow variables is derived in the subsection llll Dl In section 
IIVI we derive analytical expressions for the P.D.F., the anomalous diffusion exponents and the generalized diffusion 
constants; we also extend our calculations to a weakly dissipative system and compare our results with numerical 
simulations. Section^is devoted to concluding remarks. Technical mathematical details are given in the appendices. 
In particular, in appendix IdI we discuss the 'Best Fokker-Planck Equation' (B.F.P.E.) approximation 0,lD|, based 
on partial resummation to all orders in the correlation time r of a perturbative expansion for small noise amplitude. 



II. REVIEW OF EARLIER RESULTS 



Consider a nonlinear oscillator of amplitude x(t), trapped in a confining potential U(x) and subject to an additive 
noise : 

The statistical properties of the random function need not be specified at this stage. We restrict our analysis to 
the case where the potential U is an even polynomial in x for \x\ — > oo. A suitable rescaling of x allows us to write 

U with n > 2 . (2) 

2n 

As the amplitude x(t) of the oscillator grows with time, only the asymptotic behavior of U{x) for |x| — > oo is relevant 
and Eq.Q reduces to 

I^W+sC*) 2 "- 1 =£(*)■ (3) 

We emphasize that all stochastic differential equations involving white noise are interpreted in the Stratonovich sense, 
that allows the use of ordinary differential calculus for the change of variables. We review here some previously derived 
results that are required for the present work. 



A. Energy-angle coordinates 

A key feature of Eq.Q is that its underlying deterministic system (obtained by setting £ = 0), is integrable. The 
associated energy and the angle variable are given by [a, |6| 

Z7 1 -2 , 1 2n , j. V™ f X/El/2n du , A \ 

E= 2 X+ Yn X ' and *=(W^J y=|' (4) 
where the angle <j> is defined modulo the oscillation period 4K n , with 




FIG. 1: Equipartition relations with Ornstein-Uhlenbeck noise. Equations an d lEU are integrated numerically for 23 = 1, 
t = 5, with a timestep <5£ = 10 -5 and averaged over 10 4 realizations. Fig. (a): the first equipartition ratio (E(t)} / (v(t) 2 ) is close 
to the theoretical value given in Eq. 1121 : 3/4 for n — 2; 2/3 for n = 3; 5/8 for n = 4. Fig. (b): the second equipartition 
ratio (v(t) 2 )/{x(t) 2n ) is close to 1 for n = 2,3,4. 



In terms of the energy-angle coordinates (E, </>), the original variables (a;, x) read as follows: 

x = E 1 '^ Sniff,), (6) 
x = {2n)^E 1 ' 2 SM , (7) 

where the hyperelliptic function S n is defined as 0, E3 



The presence of external noise spoils the integrability of the dynamical system © but does not preclude the use of 
(E, (p) instead of (x, x) as coordinates in phase space. Introducing the auxiliary variable f2 defined as 



ft = (2n)— £5 ; (9) 
equation J2J is written as a system of two coupled stochastic differential equations @ 

n = nSM)m, (io) 



.(2n) 

This system is rigorously equivalent to the original problem (Eq. and has been derived without any hypothesis on 
the nature of the driving force £(t) which may even be a deterministic function or may assume arbitrary statistical 
properties. This external perturbation continuously injects energy into the system. We have shown in |fj that, 
when £ is a Gaussian white noise or a dichotomous Poisson noise, the typical value of £1 grows algebraically with 
time. We have also verified numerically that the same behavior is true for an Ornstein-Uhlenbeck noise. As seen from 
Eq. Qll|). the growth of the phase <p is faster than that of f2. Assuming that in the long time limit (f> is uniformly 
distributed over a period [0, 4K n } , we obtain, on averaging Eqs. (|6I7(I over the angle variable, the following equipartition 
relations 

(E) = =±I (i 2 > , (12) 
2n 

(i 2 ) = (* 2n >. (13) 



These relations are in excellent agreement with the numerical simulations as shown in Fig. 1. 



B. The white noise case 



In this subsection, we consider £(t) to be a Gaussian white noise of zero mean value and of amplitude T>: 

W)) = o, 

mat')) = v6(t-f). (w) 

As described in an averaging procedure over the angle variable <fi allows us to derive a closed effective Langevin 
equation for the slow variable f2. Starting with the full Fokker-Planck equation for the evolution of the P.D.F. 
Pt(£l,4>) associated with the system (I1UI11H and integrating out the angular variable 4>, we obtain a phase-averaged 
Fokker-Planck equation for the marginal distribution Pt(O): 

dtP = % (%P - - dS V with V = . (15) 

2 \ n \l J n + 1 

The effective Langevin dynamics for the variable O is thus 

{l = lh + ^> (16) 

where T> is the amplitude of the effective Gaussian white noise It is readily deduced from Eq. (|16f) that the 

variable f2 has a normal diffusive behavior with time: 

The averaged Fokker-Planck equation Qlfijl is exactly solvable, leading to the following expression for the energy 
P.D.F. 

1 1 f(2n)^E\ n ir f (2n)^ e\ , s 

p >^ = tvw)zHst) (18) 

r(.) being the Euler Gamma function The asymptotic time dependence of all moments of the energy, amplitude 
or velocity is calculated analytically from Eqs. @ El El and I18|) . leading to the long time scaling behavior of 
dynamical observables as well as analytical expressions for the prefactors (generalized diffusion constants) 

(E) = ft, 

^ = iT+T ' 



^ = F7Ttl^TT^) ■ (19) 



2n 

The physical observables grow algebraically with time and the corresponding anomalous diffusion exponents depend 
only on the confining potential at infinity. These quantitative results are in excellent agreement with numerical 
simulations Q. 

C. Scaling behavior for colored noise 

We now take £(t) to be a colored Gaussian noise with non-zero correlation time r and generated from the Ornstein- 
Uhlenbeck equation 

^f=-i« < )-i#), (20) 
at t t 

where rj(t) is a Gaussian white noise of zero mean value and of amplitude T>. In the stationary limit, when t, t' 3> t, 
we have: 

(£(*)) = and mm) = ^ e"l*-*'l/r . (21) 

z r 



-- <E> = t/2 
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FIG. 2: Asymptotic behavior of the nonlinear oscillator with additive Ornstein-Uhlenbeck noise. Eqs. Q and JUJ arc 
integrated numerically for T> = 1, n = 2, with a timestep St = 10 -5 and averaged over 10 3 realizations. For r = 0.05, 0.5 and 
5.0, we plot the average (E) vs. time t. The dashed lines correspond to the analytical predictions, (E) — 2? t/2 for white noise 
(see Eq. JTHJ) and (E) = 0.533 (£>i/-r 2 ) 2/3 for colored noise with r = 5.0 (see Eq. JBTJ). 



In [fj, we derived from a self-consistent scaling Ansatz that, when the noise has a finite correlation time r, the variable 
n has a sub-diffusive behavior and grows as 




(22) 

From this equation and Eqs. IT2l and IT3|) . the following scaling relations are deduced 

E ~ 




(23) 

We observe that the exponents associated with colored noise arc different from those obtained for white noise. We 
emphasize that Eq. (|23|) gives only qualitative scaling relations that were obtained only by heuristic arguments. In 
particular, the prefactors (which are dimensionless pure numbers) are not known at this stage. In the following 
sections, we solve this problem analytically: we calculate the P.D.F. of the anharmonic oscillator subject to Ornstein- 
Uhlenbeck noise and derive explicit and complete formulae for the moments of the dynamical observables in the long 
time limit. 

The result of Eq. ill'ol) is supported by numerical simulations: the colored noise scalings are observed for large 
enough times, even for an arbitrarily small correlation time r (see Fig- EJ • The crossover between the white noise 
scalings (|19|l and the colored noise scalings l|23() occurs when the period T of the underlying deterministic oscillator 
(which is a decreasing function of its amplitude) is of the order of t, i.e., for a typical time t c ~ V^ 1 f-^/in-V) _ p or 
t « i c , the angular period of the system is large as compared to the correlation time of the noise which, therefore, 
acts as if it were white; for t 3> t c , the noise is highly correlated over a period and is smeared out, leading to a slower 
diffusion. 

Some interesting qualitative results have thus already been obtained. But, a quantitative study of the nonlinear 
oscillator with Ornstein-Uhlenbeck noise has remained out of reach. The reason is that the averaging technique, 
that applies perfectly for the white noise, leads to erroneous results when applied to the colored noise. We now explain 
why the naive averaging method fails and how to modify it to make it work. 



III. AVERAGING METHOD FOR ORNSTEIN-UHLENBECK NOISE 



A closed Fokker-Planck equation for the probability distribution function Pt(Q, <f>) on the system's phase space can 
not be derived when the noise is colored: because the noise is correlated in time, the system's dynamics acquires 
a non-Mar kovian character and the evolution equation for P t (Q,(j)), which can be derived using functional methods 
0- El EJ 0j 13 E3- involves correlation functions at different times. This hierarchy can not be closed at 
any finite order unless truncation approximations are used |l8i l2l| . However, by embedding the initial problem 
into a higher dimensional Markovian system, a rigorous description in terms of a Fokker-Planck evolution equation 
can be formulated. Indeed, the random oscillator described by Eqs. (|10lll|l together with the Ornstein-Uhlenbeck 
noise generated by Eq. I|20(l forms a three dimensional stochastic system driven by a white noise r)(t). The stochastic 
equations 1(1011 ll and l20|l define a random dynamical system in the three dimensional space (fi, (f>, £), the three variables 
being treated on an equal footing: 

Q = nS' n (<f>)£, (24) 

o ^ (25) 



i = ~-t--v(t), (26) 

T T 

where i](t) is a Gaussian white noise of amplitude T>. The Fokker-Planck equation for the joint P.D.F. Pt(fl, 4>, £) is 
then given by 



■at = - n m{ s - {<t>HPt ) " a* {{{^) " IT ) Pt ) + r^f + J^W 



(27) 



If an exact solution of Eq. 1(27(1 could be found, the P.D.F. on the original phase space (f2, <fi) will then be obtained 
by integrating the full solution over the variable £. Unfortunately, this does not seem to be feasible. However, the 
evolution of the angular variable <fi is fast as compared to that of fl irrespective of the correlation time of the noise. 
Therefore, it is desirable to average out (f> from the dynamics in order to obtain an effective Langevin equation for the 
slow variable f2, as has been done in the white noise case. 



A. Failure of straightforward averaging 

We now average the Fokker-Planck equation (|27[) over the angular variable <p assuming that for t — > oo the probability 
measure for <p is uniform over the interval [0, 4A"„]. Using the fact that the average of the derivative of any function 
with respect to (f> is zero: 



3 (...) = O, (28) 

the phase-average of the Fokker-Planck equation (|27[) governing the evolution of the marginal distribution Pt(£l, £) is 
found to be 

dP t _ i d£p t v d 2 P t 

dt t <9£ 2 r 2 <9£ 2 ' 1 ' 

this equation corresponds to the following Langevin dynamics for the couple of variables (£^,£) 

tl = 

i = -k--v(t)- (30) 
r t 

This result implies that f2 has a trivial dynamics: it is no more stochastic and is conserved. The integration over 
the angular variable averages out the noise itself and leads to conclusions that are blatantly wrong. This failure is 
due to the destructive interference between the fast angular variable and the noise: a simplistic adiabatic elimination 
of <p eliminates the noise too. The failure of this naive averaging procedure can also be understood from a direct 
analysis of Eqs. (|24l 1251 and 126(1 : the functions S n ((f) and S' n ((j)) are periodic in <f> with period AK n and have a zero 
mean value. Thus, if we average Eqs. 1241 and 12*6^1 over (j> and discard the correlations between 4> and £, the erroneous 
dynamics (|30|l is obtained. A correct averaging scheme must take into account the fact that the colored noise £ is 
not constant but varies during the period T of the underlying deterministic oscillator (the angle <fi covers the interval 
[0,4_R"„] during T); when T <C r, the variation of £ can be estimated from its time derivative that contains a white 
noise term (see Eq. (|2(JH 'I which will survive the averaging procedure. 



B. A heuristic calculation 



The preceeding discussion suggests that we must transform the system l2~5l and |2U so as to make the time 
derivative of £ appear explicitly in the evolution equation 1|24|) of the slow variable 0. Before performing this transfor- 
mation in a systematic manner, we explain the heuristics of the method. Taking into account the fact that Q grows 
to infinity with time, let us simplify Eq. (|25(l to 

n 



■(2n) 

and use this relation to transform Eq. I|24[l as follows 



Introducing a new variable Y defined as 



Y = ^-^r, (33) 



we rewrite Eq. I|32|) as 



2n-l "^'^ di 
Using Ornstein-Uhlenbeck's equation l|26|l . Eq. i|34|) becomes 



(2nY 



1 —Y = SMM=^4^-S n (< t >)i. (34) 



V ■«)£) = ^^+«S n (^. (35) 



di v 2n-l " vr ^ t ,ivrj T 

In the long time limit, the term 5 n (0)£ which remains of finite order becomes negligible as compared to Y that grows 
as a power-law with time. Besides, on the right hand side of Eq. (|35J) we also neglect the colored noise £ with respect 
to the white noise r\. We thus conclude that Y has a normal diffusive behavior and scales as 

<F 2 )(x^. (36) 

This scaling result is equivalent to Eq. (|22|l and implies Eq. ()23|) . The variable Y undergoes an effective diffusive 
dynamics and plays a role similar to that of Q in the white noise case. In order to obtain a quantitative agreement 
between theory and numerics, we must set this heuristic calculation on a sound basis. In particular, some terms that 
were overlooked in the sketchy discussion above contribute to the long time behavior of the P.D.F. in phase space and 
must be taken into account. 



C. Coordinate transformations 



A correct averaging procedure is performed by defining recursively new sets of dynamical variables on the global 
three dimensional phase space (fi, <f>, £). These new variables will embody the correlations between <f> and £ and 
provide a book-keeping device that transforms the heuristic calculation described above into a systematic perturbative 
expansion. We show in the Appendix B that a suitable set of variables is given by (Z, </>,£), where Z is defined as 



rY^r 

the variable Y was defined in Eq. (|33|l and the function C n (4>) satisfies 



z(Y M = y - ,2„ - Dis.m + (gg^M + « . (3r , 



dC n (cj>) 



4A'„ 



C' n {4>) = ^ = $n(4>) and C n {<t>) = / C n {<t>)d<t> = , (38) 



the overline indicates an average over the angular period- 
Equation 1)37)1 defines a change of dynamical variables between the set (il, <fi, £) and the set (Z, <j>, £)■ The system (1241 
1251 and l2"F)|l must now be expressed in terms of the new variables and we obtain 

-Z = J Z (Z, 0,0+^(^,^0— , (39) 



2n-l — --v-.ri-w r 

4> = Mz,<t>,0, (40) 

£ = --€--r?(*). (41) 
r r 

The current and diffusion functions that appear in Eqs. 1)391 and l4T)|l are given by 

(2n)^feTC n (^)e , 2Sl{4>) + (n-\)S' n {4>)C n (cj ) ) t . 2 
Jz{Z,<P,0 = — 1" ~ £ . (42) 

T>(7A?\ c , (2n)^TTC»(0) , Sl(#)j 

Vz{Z,<p,0 = o n (0) H — 1 = — , (43) 

tZ^t: Z 

MZ,(p,0 = I ; - , - g. (44) 

V (2n)^ ' S2(Z,0,O 

In Appendix B, we show that the stochastic dynamical system 1)391 1401 and 141(1 is obtained from the original equa- 
tions ((241 1251 and 126)1 by retaining all contributions up the order 0(Z^ 1 ). The crucial difference between the two 
systems appears when Eq. 1)24)1 is compared with Eq. 1)39)1 : the evolution of Z contains an explicit white noise con- 
tribution that embodies the variations of the colored noise £ during a time scale shorter than the correlation time r. 
This white noise term will survive the averaging process and will allow us to derive a non-trivial Langevin dynamics 
for Z. 



D. Averaged Langevin system 



As for the white noise problem, we shall eliminate the fast angular variable in the Fokker-Planck equation for the 
P.D.F. H t (Z, </>,£) associated with Eqs. 14*01 and EFT) . This Fokker-Planck equation is given by 

g_ _ (2 „_ 1) ^ (A n,)-| wn ,, + I& + 

- V&mW - (2» - UsggPW) - P» - + ^} . (45, 

Integrating out the fast variable 4>, we obtain the following averaged Fokker-Planck equation that describes the 
evolution of the marginal P.D.F. Il t (Z,£) (the details of the calculations are explained in Appendix C) 

— = -(2n - 1) M _ (^^^n t J + - _ + (46) 

V I to ^ d2tlt to u 92 (&*\ (o u d ( Z dfL t \ d 2 fl t \ 

2^\ {2n - l) ^- {2n - \—)~ ( - l ^dz\-z^) + ^e)' (47) 

where the parameter /i, calculated in Appendix A, is given by 

— T , ,1 l2nlM 2n , „ 

M = ^ = (2n)w V 7 V 7 ^ 

r(^)r(^) 

Let us now consider the following Langevin system for the variables Z and £: 

Z = (2n - 1)M (3 ~^ 2 + (2n - l)^(t) + (2n - > ( 49 ) 

tZ t tz 

i = -U~-V(t), (50) 

T T 



where fj(t) and rj(t) are two independent Gaussian white noises of amplitude T>. Writing the Fokker-Planck equation 
related to this stochastic system I5()[) we retrieve the averaged Fokker-Planck equation (|47(1 up to terms of the 
order 0(Z~ 2 ) (recall that the stochastic Eqs. jSU E01 and ETJ are equivalent to Eqs. ((211123 and EEJ at order ©(Z" 1 )). 
We have thus found the effective stochastic Langevin equation l|49ll for the slow variable Z coupled with the Ornstein- 
Uhlcnbeck noise £. This dynamics is however not exactly solvable and we need one last transformation before deriving 
explicit analytical results. In terms of the new variable 



e 

Zj\ = /j + (zn — 

Eq. ll4*§|l becomes 



Z x = Z + (2n - , (51) 



Z X = Z + (2n - = (2n - 1) ^ ^ + (2n - l)^fj(t) , (52) 

where we have neglected all terms of order strictly less than 0(Z^ 1 ). Recalling that £ 2 has a finite mean value equal 
to T>/2t, we rewrite the Langevin equation (|52|l for Z% as follows 

Z X = (2n - 1)(2 - n)gi- + (2n - l)^(f) + (2n - 1)(2 - - (O) • (53) 

This equation contains two independent noise sources: a white noise contribution, f)(t), and a non-Gaussian colored 
noise term, (£ 2 — (£ 2 )), of zero mean and of finite variance T> 2 /2t 2 . This colored noise is multiplied by a factor 
proportional to l/z*i and therefore its effect becomes negligible as compared to that of fj(t) when t — > oo. We thus 
simplify Eq. JHSJ to 

Z X = (2n - 1)(2 - n)^±- + (2n - l)^f}(t) . (54) 
2t z Z\ t 

We have thus obtained an effective white noise Langevin dynamics for the variable Z\ . This equation plays the same 
role for the nonlinear oscillator subject to a colored Ornstein-Uhlenbeck noise as that played by Eq. Ijlfijl for white 
noise. The mathematical structure of the two equations and (|54(l is the same and they differ only by the constant 
prefactors that embody the specific parameters of each problem. 



IV. ANALYTICAL RESULTS IN THE LONG TIME LIMIT 



In this section, we solve the effective dynamics for the variable Z\ and deduce statistical results for the long time 
behavior of the energy, the amplitude and the velocity of the nonlinear oscillator driven by an Ornstein-Uhlenbeck 
process. We shall not only retrieve the scalings predicted in Eq. (j23|l but also derive explicit formulae for the prefactors 
and for the skewness and the flatness of the energy P.D.F. All these analytical results compare favorably with numerical 
simulations. Finally, we show that our analysis can be extended to include a weak dissipation term in the dynamics. 



A. Calculation of the P.D.F. 



The Fokker-Planck equation associated with the effective Langevin dynamics (|54|l of Z\ is exactly solvable (using 
the method of and we obtain the following expression for the asymptotic P.D.F. 

Pt[Zl) = r( "+i ^ (^n-in-Dt) ^ ^ 6XP ( " 2(2n-f)^2?t) ' (55) 

\2(2n-l)J 

Expressing the energy fiasa function of Z\ with the help of Eqs. 1331 PT71 and l5T | . we obtain at the leading order 




and thus dE 




dZ 1 . 



(56) 



From Eq. I|55[) . we deduce the asymptotic expression of the P.D.F. of the energy 

V ; r (ji±L) V(2n- VfyLVt) V (2n-l) 2 M X>i/ V 7 



4n-2 



The statistical behavior of the position, velocity and energy of the system in the long time limit can now be derived 
by using the P.D.F. (JSJJ and Eqs. (JOIandQ: 



1 \^4n-2 ^ 

2 "' (59) 



(i 2 ) 



n + 1 

' n+3 



^ = Fkrf( - 2n ^ ) • (60) 

1 \^4n-2 ^ 



In the last equality we used /i = <S 2 , according to Eq. I|48|). In particular, we find 
For n = 2, 



/T)t\ 2 /3 /T)/\2/3 /"D/\l/3 

(£0 - 0.533 (^) , (x 2 > = 0.711 (^) , (x 2 ) = 0.587 (^) . (61) 



For n = 3, 



For n — 4, 



, , /Vt\3/5 /Vt\3/ r o /X>i\i/5 

(£) - 0.474 (—J , (x 2 ) = 0.711 (-3-) , (x 2 ) = 0.535 (^J . (62) 



(E) = 0.436 (^) , (x 2 > = 0.698 (^) , (x 2 ) = 0.500 (^) . (63) 
From Eq. lf53)l. the skewness and the flatness of the energy P.D.F. can also be calculated 



,„„, r 7n+l \ r / n+1 

(£ 3 ) 1 (^477=2 J 1 ^47{^2 



= TT^H = — .3 ' , (64) 

r 



5w+l 

4n.-2 



(£ 2 )3/2 

, W 9?. -1 \ r / n t I 



/ 5n±l\ 
\^4n-2 y 



The skewness is approximatively equal to 1.95 for n = 2,3 and 4 and its limiting value for n — ► oo is given by 
"rflWa = r/lyz — ^-03 ■ • • ^ ne na t ness is approximatively equal to 4.7 for n = 2, 3 and 4 and its limiting value for 
n — > oo is equal to 5. 

The analytical results of Eqs. H61I65|I are compared with numerical simulations in Fig. |3| We observe that the 
agreement is remarkable in the long time limit. 

B. Extension to the case of a weakly dissipative system 
We now introduce a linear friction term in the system with dissipation rate 7. Equation l|55(l then becomes 

^ x{t)+7 A. xit)+x(t fn-l =m _ (66) 




FIG. 3: Asymptotic behavior of the nonlinear oscillator with additive Ornstein-Uhlenbeck noise. Eqs. © and 

arc 

integrated numerically for T> — 1, r = 5, with a timestep St — 10 -5 and averaged over 10 4 realizations. For n = 2,3, 4, we plot: 
(a) the average (E) and the ratio (E) / (T> t / y 2 )™/( 2 ™ -1 ) (inset); (6) the skewness and flatness factors of E vs. time t (c) the ratio 
(x 2 )/(Vt/r 2 ) 1/{2n - 1) . The dotted lines in the fi gures correspond to the analytical predictions, Eqs. l|ri()|l . I|64|l and 1)650 . 



In terms of the (il, 0) variables defined in Eqs. (JBJ and |5J) , Eq. I|66l) becomes 

it = - nl ^" { f s> + nS'M)m, (67) 

(2n)« K (2n)^y \t 

We want to average out the angular variable from this system; this is possible only if the evolution of 4> is fast as 
compared to that of the energy : this condition is satisfied when the dissipation rate is vanishingly small so that it 
has a negligible effect during a period of the underlying deterministic oscillator. We proceed as in Section II I II and 
obtain after adiabatic elimination of the angle the following white noise Langevin equation for the effective variable 
Zi [defined in Eq. 

Zi = - {2n - l ^ Zl + (2n 1)(2 - n)£± + (2n l)^(t) . (69) 
n + 1 2r z Zi t 

To derive this equation we have assumed that 7T <C 1, and have used the identity S' n 2 = ^jjrr (derived in Appendix 
lA"j) . When t — > 00, the P.D.F. of Z\ tends to a well defined stationary limit P s tat(^i) which can be explicitly calculated 
by solving the stationary Fokker-Planck equation associated with Eq. I|H9|I. Reverting to the energy variable E, we 
obtain the stationary distribution function of E 

4n 2 "/T 2 \ i£=2 i-n ( 4n 2 "/T 2 E~"~ 



^ = ^t^l (2n-l)(n + W " £? " reXp i- (2n-i)(n + W - (70) 
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FIG. 4: The weakly dissipative case. Eq. H66t with Ornstein-Uhlenbeck noise is integrated numerically for 7 = 0.02 with a 
timestep St = 10 -5 , and averaged over 10 3 realizations. For n = 2,3,4, we plot: (a) the average {E) (main plot) and the 
ratio {E)/(V/-yr 2 ) n/{ - 2n - 1) (inset) vs. V/jt 2 ; (6) (x 2 ) (main plot) and the ratio \x 2 ) / (V / 1 r 2 ) 1/{2l% - 1) (inset) vs. V/-yr 2 . 
The parameters r and T> vary within the ranges 2 < t < 10, 10 1 < V < 10 6 . The dashed lines in the fi gures correspond to 
the analytical predictions, Eqs. I|71I473|I . Averages are evaluated beyond the dissipative and colored timescales (t > I/7 and 



t > D r 



-2u/(n-l) 



respectively). 



We have thus derived, in the limit of vanishingly small dissipation, an analytical expression for the stationary proba- 
bility distribution function of the energy of a non-linear oscillator in presence of an additive Ornstein-Uhlenbeck noise. 
In contrast to the case of white noise, this stationary P.D.F. is not the canonical Gibbs-Boltzmann distribution. 

The statistical behavior of the position, velocity and energy of the system in the long time limit can now be derived 
by using the P.D.F. (|7UJ) and Eqs. (Eland©: 



3n+l \ 



( Ji±1 -) 

\4n-2 J 



An 2 "/ 



(i 2 ) = ^-(E), (72) 
n + 1 



(x 2 ) 



, p M F (i£l) / (2n-l)(n+lV 2 »P x^T 
M } ~r(n±i.\\ ) ■ (?3) 

1 \in-2 ) 



We observe that Eqs. J7U|to[73J) can be obtained from Eqs. l|57l to R)UI) by making the substitution ^ = 2 " ^ 7. 

These analytical formulae compare favorably with numerical simulations (Fig. ^| for parameter values such that 
both 7 and the product jt are sufficiently small (e.g., 7 = 0.02, 2 < r < 10). As expected, Eqs. (|71l to l73|) become 
less and less accurate as the dissipation coefficient 7 becomes larger (however, the anomalous diffusion exponent 
n/(2n — 1) in Q71JI remains roughly valid at least up to 7 ~ 1). For smaller values of r, care should be taken to ensure 
that measurements are indeed performed beyond the crossover time t c where white-noise-like dynamical behavior is 
replaced by the time-asymptotic regime characteristic of colored noise (t c ~ V- 1 T -2"/(»-i) i s a rapidly deer easing 
function of r). 

Lastly, our result for the stationary P.D.F. in the weakly dissipative regime, (Eq. I|70|l ). compares favorably with 
the P.D.F. of a linear oscillator with Ornstein-Uhlenbeck noise |22| 

'-<*■*> - ^TO 1 - { - £(t + s >) • <»> 

In the limit r> 1 and 7T <C 1, this expression reduces to 

, . 27T 2 / 2"/t 2 E\ , s 

Pstat (E) = -L- exp ( - -J^-J , (75) 



which agrees with Eq. I|70[) for n = 1. 



V. 



CONCLUSION 



We have studied the long time asymptotic behavior of a non-linear oscillator subject to a Ornstein-Uhlenbeck 
noise. Because of the nonlinearity, the intrinsic time scale of the system becomes, in the long time limit, smaller than 
the correlation time of the noise. In this regime, the noise is intrinsically colored and the behavior of the system is 
radically different from that obtained for the white noise. The difference between white and colored noise driving 
can be quantitatively probed by calculating growth exponents that characterize the long time behavior of physical 
observables such as the energy, mean-square amplitude or velocity of the oscillator. We emphasize that the colored 
noise scalings are always satisfied in the long time limit however small be the correlation time r of the noise. However, 
when r is small, a crossover from the white noise scaling regime at short times to colored noise scalings at long times 
is observed. The crossover time t c diverges to infinity as r — * 0. In other words, the limits r — > and t — > oo (t being 
the observation time) do not commute. 

In this work, our analysis is based on the exact Fokker-Planck equation obtained by embedding the initial problem 
into a higher dimensional Markovian system. In the long time limit, a dynamical separation of time-scales occurs: the 
angular variable becomes a fast variable whereas the energy evolves on a slower time scale. It is therefore tempting 
to eliminate the angle and derive an effective stochastic dynamics for the slow variable. Straightforward averaging 
leads to excellent results for the white noise case but it fails for colored noise. Indeed, at the lowest order the noise is 
practically constant during an angular period: thus, when the angular degree of freedom is integrated out, the noise 
itself is eliminated. A more sophisticated averaging procedure is therefore required that embodies the subdominant 
variations of the noise during the short time scale. By making recursive coordinate transformations on this higher 
dimensional system, we have managed to integrate out the fast variable while retaining the effect of the noise. We have 
finally derived an effective colored Langevin dynamics for the slow energy variable. At leading order, this effective 
dynamics is solved leading to an explicit formula for the probability distribution function of the energy in the long 
time limit. From this P.D.F. we have derived the moments of the position and the velocity of the oscillator. Our 
analytical results are compared with numerical simulations and the agreement is excellent. We may thus conjecture 
(without proof) that the analytical expression (|57|l for the P.D.F. of the energy, derived by our averaging scheme, is 
an exact asymptotic result. 

We have extended our analysis to the case where the system is weakly dissipative; again the agreement between 
analytical and numerical results is very good. When the dissipation rate becomes higher, the agreement becomes 
less and less accurate; indeed, the angular variable can no more be considered to be fast as compared to the action 
variable and therefore it is not possible to use an adiabatic elimination procedure. However, in the regime of very 
large dissipation rate, the inertial term can be adiabatically eliminated [23L l24l l25| and the problem is reduced to a 
first order Langevin equation with colored noise; many specific approximation schemes have been developed to study 
this case |2f|E3- The mathematical study of the intermediate dissipation regime remains to be done. Moreover, the 
extension of these averaging techniques to systems with spatio-temporal noise [28[l29j| is a challenging open problem. 



In this Appendix, we derive some mathematical relations about the hyperelliptic function S n (<p) that have been 
used in the text. In Eq. JSJ), we defined the hyperelliptic function S n as 



APPENDIX A: SOME USEFUL MATHEMATICAL IDENTITIES 




(Al) 



From this equation, we derive the following relation between S n and its derivative S' n 




(A2) 



These hyperelliptic functions reduce to circular and elliptic functions for n — 1 and n — 2, respectively: 



<Si(</>) = \/2sin0 



and 




(A3) 




and 



(A4) 



where sn, cn and dn are the classical Jacobi elliptic functions of modulus 1 / y/2 0, . 

The functions S n and S' n are periodic functions of period K n , defined in Eq. J^J. Their average over a period 
vanishes. We also need to calculate the averages of the square of these functions. Using Eq. lAljl and Eq. (JJJ, we 
derive Eq. (gSJ 

K r 1 u 2 du rf_L^p/'ii+i 

/ovfe r(^)r(^ 



To obtain this identity, we made the change of variables u = S n (<p) and then r = u 2n , and finally used the Eulerian 
integral of the first kind ^3] 

'■•(H'^ r y";" w 

1 (a + p + 2) 

Similarly, from Eqs. l|A2jl and and again making the change of variable u = S n {4>), we deduce that 



S'n = 1 . (A7) 

Using the following identity (that can be proved by integrating l.yfl — u 2n du by parts) 

/ du y/ 1 — u 2 " — n du — = — n / du \/l — w 2 ™ + n 

Jo Jo VI - u 2 ™ 7 7 



d« 



(A8) 



we conclude from Eq. I|A7|I that 



Si? = (A9) 
n + 1 

APPENDIX B: TRANSFORMATION OF THE ORIGINAL LANGEVIN EQUATIONS 

In this Appendix, we perform a change of variables in the original Langevin equations and derive the stochastic 
dynamical system l|39l EH1 an d I41f> starting from the original equations 1241 12*51 and I26|) . In terms of the variable Y 
defined in Eq. (|33|l . Eq. 124|l can be written as 

y 0, a , r ^ 



2n-l ri\(2n)sr/ V(2 n )2Tr 
Similarly, substituting the variable Y in Eq. (|25|l we obtain 

<«> 

From Eqs. (|Bll and lB2|) . we deduce 

y smk (m] 

2n - 1 1 - Ml ■ 1 ; 

1 y 

At long times, we know that Y — > oo; thus we make a perturbative expansion of Eq. I)B3J| as follows 

-Sn(<P)<t>Z+ ^7 + ^9 + ••• ( B4 ) 



2n-l »^,rs y y2 

In the following calculations we retain terms up to the order OiY^ 1 ) only, because Y — * oo in mean. In principle, 
random fluctuations may allow for small values of Y for arbitrary large times. However, in our numerical calculations, 
we never observed such fluctuations (of course, numerical simulations are carried out over a finite duration of time 



and can not be used as a proof that Y — > oo pathwise). Thus, the effective dynamics that we shall derive will be valid 
for large values of Y and the P.D.F. that we shall obtain may not describe correctly the distribution of small values of 
Y. Nevertheless, since the relative probability of observing small values of Y tends to as time grows, the formulae 
for the long time behavior of the moments will be asymptotically exact when t — > oo. 

Because we retain only contributions up to the order 0(Y~ 1 ), we can neglect the last term of Eq. {B4j: indeed, from 

Eq. I|B2(1 we observe that this term is of order o{^^j = 0{Y" 2n ~ 1 ) i which is less than 0(Y _1 ) . We now integrate by 

parts the two relevant terms on the right hand side of Eq. (|B4Jl in order to make the white noise contributions appear 
explicitly: 

2^1 " dt \ SnW ^ + ~^Y~ ) + r + tY + Y^ + ■ ' ' (B5) 

The last term in this equation is of the order O^y^j = 0(^2 J and is negligible with respect to ©(Y -1 ) . We thus 
obtain up to the desired order 

dt\ 2^~l SnW * -^Y-J = r + tY ■ (B6) 

On the right hand side (r.h.s) of this equation, the white noise terms (proportional to rj) and the term <S> 2 (0)£ 2 / VY will 
survive the averaging over the angular variable <$> whereas the term <S„(0)£/t will be eliminated because S n (<p) = 0. 
However, this term contains important correlations between £ and <j> of the order 0(Y^ 1 ) and in order to retain these 
correlations, we transform this term as follows 



V(2n)3n/ V 



Y 



= — 1 5^~2 1- • • • (,B7j 

Integrating by parts the first term on the r.h.s. of Eq. I)B7|) and using the function C n ((f>) defined in Eq. I|H8|I. we 
obtain 

(2n)^^iy«S n (^)^ _ d /(2/.)sfelC n (^)^ , (2^^,^) (£ + r,) , (2n)5i^J(n - !)£„(<£) £ Y 



-rY^r-i dt\ tY ^t I T 2 Y 2^r (2n-l)rY^=r 

d / (2n)^r>C» (</,)£ x | (2n)^WC ra (0)(£ + r ? ) | (n - 1) S> n (<f>)C n (0) £ 2 
diV T yir^T / T 2yi^i rY 

where we have used Eq. (jBlfl to derive the last equality. Similarly, the second term on the r.h.s. of Eq. I|B7|) is 
rewritten as 

(2n)^s^e = me l = me + 0(r _ 9) t m 

T y5^i tY f Y \ss=t tY 



where we have used Eq. l|B2|l . If we substitute Eqs. l)B8l and tB9jl in Eq. I)B7|1 and then substitute the result in Eq. (|B6(I . 

we derive the following equation for Y, valid up to the order C(Y _1 ) 

d f Y (2n)^T)C Tl (^)£ 5|(0K 2 _ l 

tH q 7 ~<5n 4 ^ ^ f = 

dr L 2n — 1 T y— 2r J 



( 2n )^(^T)C»(0)e | 2S 2 (0) +(n-l)5;(0)C n W 
r 2y^r tY 



T 2 Y 2T^T 

2, w a W + fi!ggg^ + glMi}, (B10) 



2 



The left hand side of Eq. (|B10j) is identical to 2^rT: wit h ^ defined in Eq. p7|. On the right hand side of Eq. l|B10|) . 
we must express Y as a function of Z. Although the relation Ij37|l between Z and Y can not be inverted by a closed 
formula, Y can be calculated as a function of Z up to the order 0(Z~ 1 ) included. This leads to Eq. I|39|) and to the 
formulae (|42l and l43fl for the current and diffusion functions. 



APPENDIX C: DERIVATION OF THE AVERAGED FOKKER-PLANCK EQUATION 



In this Appendix, we derive the averaged Fokker-Planck equation l|47|) by eliminating the fast angular variable from 
Eq. H45J) . We assume that the P.D.F. II t (Z, <fi, £) becomes independent of <f> when t — > oo, i.e., the probability measure 
becomes uniform with respect to the angular variable in the long time limit. We integrate Eq. H45J) with respect to 
the fast angular variable 4> and examine the behavior of each term of Eq. I|45[) . The average of the current Jz(Z, 4>,£) 
is given by: 



Jz = — — 1 £ • (.^1) 



Integrating by parts and using Eq. H38|) . we obtain 



s>M)Cn^) = -s n ^)C'M) = -s*(4>) 

Using Eq. (|48(l . we deduce that 



Therefore, we have 



-z = ^»e. (C2) 



Similarly, the average of the diffusion term T>z{Z, is given by: 



V z =S r M + — ^r^ + ^f^ = ^, (C4) 

tZ~ Z / 

where we have used S n = C n = 0. The angular average of the expressions (T^z^t) and -§^D z^gjh that appear in 
Eq. i(4"5|) are readily deduced by using Eq. (|C4|) : 



Jk ( ^ nt) = (C5) 



az p ^ = "az^zarj" ( } 

Recalling that ^ scales as Z -1 , the leading order in the expression -§^0 'z -§z^P 'z^lt) scales as Z~ 2 . Therefore, 
retaining only the leading term in the average of this expression, we obtain 



Finally, using the fact that the integral of any term of the type c^(. . .) over a period of cj) vanishes, we deduce 



^(J n t ) = o. (C8) 

We have thus calculated the angular average of each term that appears on the right hand side of Eq. (|45|l [see 
Eqs. (jnHEnHniinilandinHl]- This concludes the derivation of Eq. l|t7|. 



APPENDIX D: THE 'BEST FOKKER-PLANCK EQUATION' APPROXIMATION 

In this appendix, we shall discuss an approximation called the 'Best Fokker-Planck Equation' (B.F.P.E.) which is 
based on a partial resummation of the perturbative small r expansion 

. This B.F.P.E. has been rightly criticized 



|30| because it can lead to quantitatively incorrect results (e.g., in the study of the activation rate of overdamped 
bistable oscillators driven by colored noise). Nevertheless, we shall show that, for the system under study, the B.F.P.E. 
leads to the correct colored noise scaling behavior; it can also describe the crossover between white and colored noise 
regimes as we showed in [TT|| . Because the B.F.P.E. can not be expected to yield exact quantitative results, our 
discussion will be mostly qualitative. 

Because random dynamical systems with colored noise are non-Markovian, there exists no closed Fokker-Planck 
equation that describes the dynamics of the P.D.F. in phase-space. The Markovian character can be retrieved by 
embedding the system into a higher dimensional process, as has been done in the present work (see Section . Nev- 
ertheless, many approximation schemes have been developed to derive effective Fokker-Planc k eq uations for systems 
subject to colored noise such as short correlation time expansions 0,0, E3, functional methods [lPlEol ]. unified colored 
noise approximation (U.C.N.A.)|2a|. or decoupling Ansatz |2l| (for a review, see 0,0, We have compared 

some of these approximations with our averaging technique and our numerical results. We have observed that no short 
correlation time expansion truncated at any finite order could yield the correct colored noise scalings of Eqs. 123fl . 
Because, in the long time limit, the correlation time r of the noise is not the shortest time scale of the system, any 
finite perturbative expansion with respect to r is bound to fail. An expansion that embodies terms of all orders is 
needed: an approximate resummation of this type is provided by the B.F.P.E. 

The B.F.P.E. is constructed from an expansion of the stochastic Liouville equation associated with the system i|10l 
and lll(l . These equations are rather involved due to the presence of hyperelliptic functions and because of the indirect 
coupling between £1 and <f> through the noise term £. We can radically simplify these equations without altering their 
qualitative properties by writing 

9. = K COS (f)t;(t) , (Dl) 

^(A)^- (D2) 
\ (2n) 2™ / 

We have neglected the l/f2 term with respect to f^™ -1 )/™ in the equation for <fi. Moreover, the hyperelliptic function 
S' n ((j)) has been replaced by the ordinary circular function Kcosd. The constant k will be adjusted later so that the 
white noise limit of the system IjDll and |D2 |) coincides with that of the original system IjlUI and fTTf) . 

The second-order cumulant expansion [lj of the stochastic Liouville equation corresponding to Eqs. (|10l and lTT|) is 
an approximate evolution equation for the P.D.F. P t in the small noise limit. It is given by 



m 
at 



L P t + / dx(U(t) exp(L x)Li(i - x) exp(-L x))Pi , (D3) 
Jo 



where the differential operators are defined as 

L ° p ' = -^(((i^^'). (D4 » 

U(t)P t = - ^ (f(t)/s cos <f>P t ). (D5) 
Following the procedure of ^2] > we evaluate the right hand side of Eq. (|D3J| by applying the following operator formula 



exp(A)B cxp(-A) = B + [A, B] + I [A, [A, B}} + I [A, [A, [A, B}]] + ... , (D6) 

with A = Lo and B = ~Li(t — x); we must therefore calculate iterated commutators of the operators Lo and Li. By 
induction, it can be shown that the n-th commutator is of the form 

T k = [Lo, [. . . , [L , [Lo, Li(t - x)]] . . .] - - x) f Jrff{ k) («, 4>) + ^-H^ (O, </>)) , (D7) 



dQ 1 ' ^' d<j) 



where the functions and flj* sa tisfy the recursion relations 



H 2 



{2n)z 

n—1 a--n ( H \ dHi k ~^ 
r #i - r — h — . (D9) 



and the initial conditions — — cos</>, — 0. The solution of these recursions is readily guessed once the first 
few terms have been calculated and we obtain 



H{ k) = (-iMt^t) " cos^ + fcf), (D10) 

V (2n) 2n / i 

Hr (.^(^(^^^I,. (DU) 

Using the expressions IjDIOf) and 1D11|) for and , respectively, and the autocorrelation function 121() of the 
Ornstcin-Uhlenbeck noise, we calculate the integral on the right hand side of Eq. I|D3|I 



/ da;(Li (t) exp(Lox)Li (t — x) exp(— Lox)) 
J o 



2r V 5 ^ 2 



where we have defined 



Hr{n,4>,t) = ^ j -2^ 2#> and W 2 (n, ^ t) = £ ^L! (D13) 

fc=0 ' fc=0 

In the limit t — > oo, the integral JJ^ T da;a; fe e _;l: converges to fc!, and the series defining Hi and can be calculated 
in a closed form : 

^ ,„ , \ cos0 + Tw sin <f> ,s 

Hi ft, 0,00) = t , D14 

1 + r^w z 



2 n— 1 51 »> (1 — t 2 cj 2 ) cos0 + 2toj sin< 
n (2n)^ (l + r 2 w 2 ) 2 



H 2 (ft,0,oo) - ^~ ~ - zr^h^ -• (D15) 

where we have defined 



.(2n) 

Substituting these expressions in Eq. I)D12|) . we derive the B.F.P.E. associated with the system (|Dll and lD2|l 

dP t d / / 51 v n=i \ <9 2 / cos 2 + rw sin 4> cos \ 

~dt^~ 90^(2^-' " V + ~9Q2" ^ 1 + T 2 W 2 *J 

Vkt cos 4> n — 1 9 2 / Q~n (1 — t 2 w 2 ) cos + 2tlo sin 



n 95190 \(2n)^ (1 + t 2 cj 2 ) 2 



(D17) 



Integrating out the fast variable cf>, we obtain an averaged B.F.P.E. for Pt(51), the probability distribution of the slow 
variable, 

dPt_Vnd, 1 m\ 

dt 4 951 V 1 + r 2 w 2 an J' y ' 

For t = 0, this equation leads to a scaling behavior identical to the averaged white noise Fokker-Planck equation i|15|) 
and the average of 51 2 agrees with Eq. i|17fl if we choose 

k = (2n)^ . (D19) 
For a non-zero r, the scaling behavior predicted by the averag ed B.F.P.E. (jD18fl is indeed 51 - f/(4n-2) which 

is m 

agreement with Eq. (|22|l . For rCl, Eq. (|D18|I describes the crossover between the white noise regime at short times 
and the colored noise regime at long times 
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